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FOREWORD 

This report presents results of the fourth in a series of studies 
to analyze earth orbital maneuvers including stationkeeping and dock- 
ing. The work was performed by Lockheed's Huntsville Research & 
Engineering Center for NASA-Mar shall Space Flight Center under 
Contract NAS8-30515. The three previous studies are documented 
in References 1, 2 and 3. 

The present study was performed during the period 1 July 1971 
to 6 March 1972 in accordance with Modification 9 of Appendix F. A 
second independent study — "High Frequency Cutoff Filter Analysis," 
was completed under this same contract and is documented in Volume II. 

The NASA-MSFC technical coordinators for the study were Mr. 
E.C. Smith and Mr. M. Brooks of S&E-ASTR-SG. 
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Section 1 

INTRODUCTION AND SUMMARY 

Many space missions planned for the future will require orbiting space 
vehicles to maneuver in close proximity (less than 1000 meters) to each other. 
Mission tasks requiring such maneuvers include stationkeeping, docking, orbital 
assembly, satellite inspection and retrieval, and payload deployment. Before 
a system is designed to perform these tasks, the critical parameters must be 
defined and alternative concepts investigated. Some major variables include 
mass and inertia properties, vehicle geometry, sensor data (amount and accu- 
racy), reaction control system sophistication and manual involvement in the 
control loop. 

During this study a simulation capability was developed to evaluate the 
alternative methods and hardware proposed to accomplish the tasks described. 
The fixed base simulation is performed on a hybrid computer, and includes a 
cathode ray tube display so that the maneuvers can be visually monitored by a 
pilot located in the active vehicle. Both automatic and manual control modes 
are included and the pilot can switch between these two modes at his discretion. 
Major simulation components which were developed in previous studies are: 

(1) Equations for generating the CRT display (Ref. 1); (2) two-burn translation 
logic for automatic docking (Ref. 2); and (3) baseline automatic attitude control 
system (Ref. 3). 

A complete description of the hybrid simulation is presented in Section 2, 
and simulation results are described in Section 3. Major digital subroutines 
are described and flow charted in Appendix A, and the analog wiring diagrams 
are presented in Appendix B. 
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Section 2 

PROGRAM DESCRIPTION 
2.1 GENERAL DISCUSSION 

The hybrid computer program developed during this study simulates 
the motion of a chase vehicle as it maneuvers for stationkeeping, inspection, 
or docking relative to a target vehicle. A typical docking situation is illus- 
trated in Fig. 2-1. The target vehicle is assumed moving in a circular earth 
orbit with its attitude fixed relative to some reference axis (either the local 
vertical coordinates pictured or an inertial frame). Chaser six-degree-of- 
freedom motion is controlled with attitude and translation reaction control 
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thrusters and is disturbed by gravity gradient torques. Translation and atti- 
tude maneuvers are limited by analog scaling as follows: 

Maximum Range = 1000 meters 

Minimum Range = Precontact 

Maximum Velocity = 5.0 meters/second 

Maximum Angular Rate = 0.1 radians /second 

Five major program options are listed below. 

1. Chaser Control Mode 

• Automatic 

• Manual 

• Real time switching between automatic 

and manual control at the pilot's discretion. 

2. Target Attitude 

• Local vertical hold 

• Inertial hold 

3. Sensor Accuracy 

• Perfect sensors 

• Effect of sensor inaccuracies introduced 
through measureables 

4. Application of Control Forces 

• Applied along body axes 

• Applied along local vertical axes 

• Set identically equal to zero 

5. Automatic Control Torques 

• From baseline attitude control laws 

• Set identically equal to zero. 

Computer runs involving pilot interaction (manual and switching control 
modes) must be made in real time. Runs which use purely automatic con- 
trol, however, can be made at accelerated run speeds (currently 10 times 
real time) to achieve a savings in computer time. 
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Figure 2-2 is a block diagram which shows the major elements of the 
hybrid simulation. This figure also shows how the work load is distributed 
between pilot, digital computer, and analog computer. Sections 2.2 — 2.6 

further describe the program operation and capabilities. A summary of re- 
sults is presented in Section 3. 

Major digital subroutines are described and flow charted in Appendix A, 
and the analog wiring diagrams are documented in Appendix B. 

2.2 TARGET DISPLAY AND DOCKING CONDITIONS 

A representation of the target vehicle is generated on a CRT display 
scope for visual feedback to the pilot. In the present display the target is 
represented by a cylinder. A bright spot on one of the cylinder end circles 
represents a docking target and defines the front face. Clues as to the rela- 
tive target attitude are obtained from spot location, direction of the cylinder 
centerline, flattening of the end circles, and relative size of the end circles. 

The overall display size gives range information. The pilot attempts to fly 
into the cylinder front face along the cylinder centerline which is assumed to 
be the docking axis. The angular field of view displayed on the scope is a 
program input called ctp. All display elements are normalized by before 
positioning on the screen. 

Figure 2-3 illustrates the cylindrical representation of the target. 
RHOZM, RHOZP and RADIUS are program inputs. During the simulation the 
target holds a constant input attitude relative to some reference frame (inertial 
or local vertical) as shown in Fig. 2-4. Normally, during checkout runs the 
target is kept aligned with the local vertical coordinates so that the docking 
axis (-X^,) is directed toward the earth. The chase vehicle is considered 
docked, and the run terminated automatically when the range to the target 
front face is less than some input value, PM1. 

During periods of automatic control the chase vehicle translates to a 
rendezvous point located along the docking axis and defined by the input vari- 
able PM2. When PM2 is greater than PM1, docking can only be completed 
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Attitude Rates 


Fig. 2-2 - Block Diagram of Hybrid Docking Simulation 
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Fig. 2-3 - Cylindrical Representation of Docking Target 



Fig. 2-4 - Docking Geometry 
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by switching to manual control. Nominal data for the target display and 
docking conditions is shown in Table 2-1. The input value for was com- 
puted as follows so that after a perfect dock the target front face will just 
fill the CRT display scope. 

. -1 RADIUS 

“F =tan "PM1 ' 


Table 2-1 

NOMINAL TARGET AND DOCKING DATA 


Input Name 

Nominal Value 

RHOZM 

5.0 m 

RHOZP 

5.0 m 

RADIUS 

3.0 m 

PM1 

1.732 m 

PM2 

2.0 m 

ALFAF (a F ) 

1.05 rad 


2.3 CHASE VEHICLE EQUATIONS OF MOTION 

The chaser equations of motion are solved on the analog computer. 
These equations are described briefly below and the detailed wiring diagrams 
are shown in Appendix B. 

• Translation Equations — The translation equations are solved in local 
vertical coordinates and were derived by Clohessy and Wiltshire (Aerospace 
Science , September 1960). 

x(l) = Fj + 2cox(2) + 3co 2 x(l) 

x(2) = — F 2 - 2 C0x(l ) 

x(3) = ~ F, - 0) Z x(3) 
m j 

The terms x(l), x(2) and x(3) are the local vertical coordinates defined in Fig. 
2-1. The F^ are control forces, m is the chaser mass, and CO is the target 
orbital rate in radians per second. A dot above a variable indicates differ- 
entiation with respect to time. The cross -coupling terms which are all a 
function of co can be removed from the simulation with a simple switch. 
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• Rotational Equations — Body angular rates (P,Q and R) are computed 
as shown below. Note that body axes are assumed to be principal axes. 
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Control torques, TC, and gravity gradient torques, TG, are computed on the 
digital and sent to the analog separately to facilitate scaling. 

• Kinematic Equations — The chaser attitude is computed by integrating 
a set of four Euler parameters which relate body axes to inertial axes. These 
equations have no singularities so that all attitudes can be simulated. Ortho- 
gonality is assured by forcing the sum of the squares to equal one. This is 
all done continuously on the analog. 

Nominal values for chaser mass properties, control forces, and control 
torques are shown below. 


Chaser 

Parameter 

Nominal Value 

Mass 

1.264 x 10 5 kg 

I 

4.203 x 10^ kg-m 2 

XX 


I 

yy 

2.052 x 10 7 kg-m 2 

i 

2.319 x 10 7 kg-m 2 

zz 


F. 

7.340 x 10 3 N 

i 


TC 

2.712 x 10 4 N-m 

X 


TC 

5.762 x 10 4 N-m 

y 


TC 

5.084 x 10 N-m 

z 
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2.4 AUTOMATIC CONTROL 

The control laws for attitude control during automatic docking were de- 
veloped and incorporated into the simulation during a previous study (Ref. 3). 
These control laws are sub -optimal when the minimum impulse bit is one- 
third of the digital sampling interval, but they do not include the effects of 
system time delays. When the automatic and manual modes were combined 
to permit switching, a 0.14 second system time delay was introduced into the 
automatic mode. This was necessary to ensure compatability between auto- 
matic and manual control. The time delay caused rate overshoot during the 
elimination of attitude errors and also introduced limit cycle behavior while 
the chaser was in attitude hold. 

This problem was remedied without sacrificing system accuracy by 
compensating for the time delay in the control logic. Figure 2-5 shows time 
histories of the chaser angular accelerations and angular rates with and with- 
out the control law modifications. These runs were initialized with a 10- 
degree attitude error about each of the chaser axes. The cause of instability 
and the control law modifications required to eliminate it are discussed below. 

The optimum switching lines derived in Ref. 3 are valid for only a very 
restricted number of inputs. Two categories of imperfection may arise: 

1. Errors, miscalculations, and finite time delays in the com- 
puting elements required to generate the optimum switching 
line, sense the state of the system, and switch the torque at 
the proper instant. 

2. The effects of system inputs or loads other than those used to 
derive the optimum switching line. 

These imperfections will cause slightly nonoptimum switching which, in turn, 
may permit the state of the system to overtravel the switching line, reverse, 
overtravel again, and so forth. The system state then approaches equilibrium 
while executing a high frequency oscillation in the vicinity of the switching 
line. As an example, the effect of time delay is shown in Fig. 2-6. In this 
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Fig. 2- 5 - Chaser Angular Accelerations and Angular Velocities Versus Time During the 
Elimination of 10° Attitude Errors and with a 0.14 Second System Time Delay 
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case the actual switching line is a parabola of the same shape as the optimum 
(commanded) switching line but displaced from it in both the E and E direction 
From this example we see that late switching caused by system time delays 
will produce wasteful firings by the control system. 

E 



Fig. 2-6 - Phase Plane Description of How System Time Delays Affect 
the Control Logic Switching 


The equation of the parabolic trajectory in the lower half plane follow- 
ing a step input is 


• (E 


E o> 


where I and L are the moment of inertia and torque about the axis of 
interest and E Q is the initial error. In the phase plane E = dE/dt or 
dt = dE/E. Therefore, in the lower half plane, 


E 
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and we can express E and E as follows: 


E 



+ E 

o 


E 


Lt 

I 


If we use the subscript sw to denote system conditions at the switching point, 
the switching line can be expressed as: 


also 


E 


sw 


I 

2L 


(E ) 
sw 


sw 


— (t ) 2 + E 
21 ' sw 7 o 



Lt 


sw 

I 


( 1 ) 


Now the time at which the torque reversal actually takes place is delayed 
beyond the time t gw , by the time delay, T, so that the delayed switching 
line occurs at 


or 


E' 


Lt' 


sw 


sw 


= E 


sw 


LT 

I 


E' 


sw 


L(t' ) 2 
sw 

21 


+ E = E + (E' )T + ~- 
o sw sw 21 


sw 


E' + 
sw I 


E = E 
sw 


sw 


(E 1 )T + 
sw' 21 
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Substituting into Eq. (1) 


E 


sw 


(E ) T + 
' sw 



I 

2L 



rearranging, we get 


E' + 
sw 




2LT 

I 


2 


This is the trajectory of the lower delayed switching line — with apex occur- 
ring e' sw = -Lt 2 /2 and E gw = -2 Lt/I. Similarly, the trajectory of the 
upper delayed switching line is 


E 


i 

sw 


L T 
I 



2LT 

I 


2 


Its apex occurs at 


E* = -Lt 2 /2 and £' 
sw ' sw 


2LT/I . 


The original optimum switching lines were designed for zero time delay. 
In order to compensate for the effect of a constant time delay, it was neces- 
sary to advance the switching time by an amount exactly equal to the time 
delay. We now let t” w be an advanced switching time which occurs T 
seconds prior to t gw . Then the advanced switching line occurs at 


e" 

sw 


L(t" ) T _ 

sw' _ £ , Lr 

I sw I 


L(t 


sw 


sw 


21 


+ E = 
o 


sw 


(E n )t + ~- 
v sw' 21 


or 


E 


sw 


= E 


LT 


sw 


E. 


sw 


E + (E )r - 
sw v sw 7 


LT 

21 


( 2 ) 

(3) 
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Substituting into Eq. (1) 


E 


11 

sw 


+ [}KJ r 



I 

2L 




2 


rearranging we get 


tr 11 + jlZI — I I tt" 

E sw + I " 2L [_sw 


2Lr 


— i 2 


(4) 


This is the trajectory of the lower advanced switching line with apex occur- 
ring at Eg' w = -LT 2 /I and E gw = 2LT/I . Similarly the trajectory of the 
upper advanced switching line is 



I 

2L 


E + 2LT/I 
sw 


(5) 


Its apex occurs at E^' w = *Lt /i and E gw = -2Lr/l . These advanced 
switching lines are the commanded switching lines required in the presence 
of a time delay. They are shown in Fig. 2-7 . 



Fig. 2-7 - Phase Plane Representation of Advanced Switching Lines 
to Compensate for Time Delays 
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Equations (4) and (5) reduce to the optimum switching line when r = 0 
(no time delay). 


or 


e" 

sw 

' Tu [<.] 

lower 

switching line 

e" 

sw 

* -ri;[ 6 sw] 2 

upper 

switching line 

e" 

sw 

= _ JL e" |e” 

2L sw l sw 




Implementation of Eqs. (4) and (5) is considerably different from that of Eq. 
(6) and would require a major reprogramming effort. However, for most 
cases, the time delay and L/l are fairly small so that Eqs. (2) and (3) can 
be approximated by: 


A* T-»t * 


sw 


= E 


sw 


sw 


E" + E T 
sw sw 


The approximated switching lines can be written as 


E + E T S 
sw sw 


I f" e" 

2L sw swi 


which has the same form as Eq. (6). The only change made in the raw con- 
trol law was to replace E by E + E T. 


The objective of the zero velocity control law is to eliminate limit 
cycles and to provide a sufficient buffer for the analog computer noises. 

The presence of time delay can cause unnecessarily large residue velocities 
due to excessive firing. These large residue velocities often lead to fast 
drift and make the system more oscillatory. In order to prevent this 
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shortcoming, the commanded no-firing regions inside the deadzone and in 
the first and third quadrants on the E, E plane were enlarged (see Fig. 2-8). 
The enlarged areas correspond to that velocity increment, AV, which would 
result from continuous firing for a time span of sampling period plus time 
delay. This change in the zero velocity control laws ensures that the velocity 
of the system in the second and fourth quadrants is always small, thus de- 
creasing drift and making the system less oscillatory. Note that the actual 
no-fire region remains unchanged. 

Translation control for automatic docking is accomplished using the 
two -burn translation logic developed in Ref. 2. This logic attempts to com- 
plete docking in a specified (input) time using a maximum of two burn 
periods in each axis. A separate deadzone logic holds the docked position 
once it is reached. This deadzone logic acts independently in each axis 
and requires that the chaser stay in local vertical attitude hold during 
automatic docking. This system produces minimum fuel docking transla- 
tion when coriolis, gravitational and centrifugal terms are neglected in the 
translation equations. 

2.5 MANUAL CONTROL 

Chaser control during the manual mode is modeled after the Lunar 
Module Reaction Control System used in Ref. I. Translation maneuvers 
are made by commanding control forces along chaser body axes. The 
chaser is permitted to coast when translation commands are zero. The 
attitude control logic which is described below maneuvers by rate com- 
mands. System delays are simulated by delaying all commands by 0.14 
second before their effect is felt through the thrusters. 

The chaser angular rates are compared with pilot commanded rates, 
p Q c and R c> from the attitude hand controller. The resulting error 
signals, , are fed back to the control logic. The subscript i will be 
used throughout this discussion to represent the chaser body (principal) 
axes, XB, YB and ZB. 
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E 



Fig. 2-8 - Enlarged Phase Plane Region for Zero Thruster Commands 
to Compensate for System Time Delay, T. 
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There are two modes of manual attitude control, a maneuvering mode 
and an attitude hold mode. The maneuvering mode is used whenever there 
are rate commands from the attitude hand controller i.e., whenever P , Q c> 
and R_ c are not all equal to zero. In this mode only the angular rate errors 
about the three-body axes, Ek, are fed back. These three rate errors are 
tested against a rate error deadband, as shown below. 


E. 

1 



If Ek is greater than REDB. then the manual control logic calls for a negative 
torque about the i axis. If Ek is less than -REDIk a positive torque is called 
for. The REDEk are computed internally as a function of chaser inertia char 
acteristics, control torques, and digital sample time. 

REDB. = (5.) (DEL.) 

DEL. = (TSAMP)(TC./l..) 

Note that vehicle attitude is allowed to drift in this mode. 

Whenever there are no rate commands (i.e., P c = Q c = R c = 0) the atti- 
tude hold control laws are used. In this mode attitude errors about the body 
axes are fed back along with the angular rate errors. Figure 2-9 illustrates 
the phase plane for the attitude hold mode. The parameters which define the 
torquing regions and coast region are computed as shown in the following: 
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Fig. 2-9 - Phase Plane for Manual Attitude Hold Mode 
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(Outer Rate Limit)^ = (5.)(DEL^) 

(Inner Rate Limit)^ = (3.)(DEL^) 

(Rate Crossover). = DEL^ 

(Error Crossover). = 1.4 degrees 

A new holding attitude is computed every time that the simulation switches 
to manual attitude hold. When the attitude hold mode is entered, the program 
starts comparing the three angular rate errors to the outer rate limit in the 
attitude hold phase plane. As soon as the magnitude of all three E^ are less 
than their corresponding outer rate limit, the holding attitude is set equal to 
the present inertial attitude. The manual attitude hold mode will then main- 
tain this attitude within the designated deadbands until attitude maneuvers 
are commanded. 

Note that if the attitude rate commands are equal to or less than the 
outer rate limit, then there will be no attitude coasting when the attitude hand 
controller is returned to zero. If the rate commands are greater than the 
outer rate limit however, the astronaut must allow for some attitude coasting 
after he returns to attitude hold. 

2.6 SIMULATION INPUT-OUTPUT 

All program input and output has been converted from English to metric 
units. Inputs are made in three ways — potentiometer settings, card reader inputs, 
and typewriter inputs. Chaser mass and inertia data, control torque magnitudes, 
and target orbital rate are input to a digital POT- SET routine. This routine 
then computes the required potentiometer settings and automatically sets all 
the required analog potentiometers. Data read in on cards include option 
flags, control system variables, docking geometry, sensor errors, digital 
sample time, and initial chaser attitude, velocity, and angular velocity. 

Chaser initial position, run identification and desired docking time are typed 
inputs . 
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The simulation stops automatically when the chase vehicle has success- 
fully docked or crashed. Runs can also be terminated through a manually 
controlled switch or at a time specified by typed input. Upon run termination 
the digital computer prints out the case number, the initial conditions, and 
the reason for termination — successful docking, missed attempt, or run 
flagged off. Also included in the printout are all the end conditions needed 
to evaluate the docking run, which are: 

• Run time 

• Fuel consumed and number of jet firings 
— this is presented separately for auto- 
matic mode, manual mode, attitude control, 
and translation control 

• Chaser position, velocity, orientation, and 
angular rate 

• Lateral offset from docking axis 

• Lateral velocity 

• Roll misalignment 

• Axial misalignment 

• Axial angular rate. 

Some of these docking parameters are illustrated in Fig. 2-10, and the corre- 
sponding Apollo design values are listed. 

Two strip chart recorders with eight channels each are used to record 
time histories of some of the docking variables. Usually position, velocity, 
angular rates, and jet firings are recorded. 
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Parameter 

Closing Velocity 
Lateral Velocity 
Lateral Offset 
Axial Misalignment 
Axial Angular Rate 
Roll Misalignment 


Apollo Design Values 

0.030 to 0.3048 m/sec 
0. 1 524 m/sec 
0.3048 m 
10.0 deg 
+1.0 deg/sec 
+ 10.0 deg 



Lateral 
V elocity 


LM/ATM c.m. 


Fig. 2-10 - Docking Parameters Computed After Run Termination 

J 

/ 
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Section 3 

SIMULATION RESULTS 

Hybrid simulation runs were made using the nominal target and chaser 
data to verify the program and investigate certain system parameters. Pilot 
capability to perform during manual control was limited because no hand con- 
trollers were available. Manual runs did however verify program subroutines 
and permitted a qualitative evaluation of the nominal chase vehicle for manual 
control. 

In all manual runs the pilot was located at the chaser center of mass 

with his viewing direction along the chaser x-axis. Runs showed that manual 

docking within Apollo design accuracies was possible using a crude switchbox 

(six three-position switches) for attitude and translation commands. These 

runs also indicated that the nominal control forces (7340 Newtons) are suitable 

for accurate manual translation. They produce linear accelerations of 0.06 

2 

meters per second . The nominal control torques (based on two jets per axis) 

2 

produce angular accelerations of about 0.18 degree per second . These are 
undesirably low for manual maneuvers. The low accelerations require long 
attitude maneuvering times which permit error buildup in the non- maneuvering 
axes. This puts a greater strain on the pilot. Runs with increased control 
torques (four and six jets per axis) improved manual controllability but re- 
sulted in greater fuel consumption. 

Automatic docking runs were made to show the effects of initial chaser 
position, attitude, allowable docking time, and disturbance forces on fuel re- 
quirements. In all these runs initial chaser position is expressed in the local 
vertical coordinates illustrated in Fig. 2-1. 

Figure 3-1 shows the fuel required to automatically dock from 305 meters 
(1000 feet) as a function of docking time. Initial chaser range was equally divided 
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Initial Position 
x(l) = 176m 

x(2) = 176m 

x(3) = 176m 


O Cross -coupling terms omitted 
A Cross -coupling terms included 
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Fig. 3-1- Fuel Required to Automatically Dock from 305 Meters 
Distance vs Time to Dock 
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between the three local vertical axes. Chaser attitude errors and velocity 
were initially zero. Runs were made both neglecting and including the 
coriolis, centrifugal, and gravitational disturbance forces. These curves 
show the large fuel savings which result when docking time is increased (i.e., 
coasting velocity is lowered). The effect of the disturbance terms is primarily 
to lengthen docking time rather than increase fuel consumption. The fuel used 
for attitude control increased with docking time but was always less than 7% 
of the total. Figure 3-2 presents the time histories of chaser position and 
velocity for the last pair of points in Fig. 3-1. As expected, the disturbance 
effects are felt in the x(l) and x{2) axes and not in the x(3) axis. 

Table 3-1 shows how docking time and fuel requirements are affected 
by the chaser approach axis. In all six runs the chaser initial range was 
300 meters and the desired docking time was 300 seconds. 


T able 3-1 

TIME AND FUEL, TO DOCK VS APPROACH AXIS 


Initial Chaser 
Displacement Axis 

Time to Dock 
(sec) 

Fuel to Dock 
(N) 

+ x(l) 

324 

3780 

-x(l) 

322 

3110 

+ x(2) 

298 

3130 

-x(2) 

296 

3340 

+ x{3) 

295 

1670 

- x{3) 

296 

1430 


Figure 3-3 shows how docking fuel requirements increase with range when 
docking time is kept constant. This figure also shows the maximum coasting 
velocity for each run. 
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Although the results shown above include the fuel required to maintain 
attitude hold, the runs were all initialized with zero attitude errors and rate 
errors. Figures 3-4 and 3-5 show the fuel required to eliminate attitude 
errors and attitude rate errors, respectively. 

Simulation runs were made using the switching option (real time switching 
between automatic and manual control). One sequence which showed promise 
was to maneuver automatically to a rendezvous point near the docked position 
and then switch to manual control for final alignment and translation. This 
method gives low fuel consumption and good accuracy. 
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Fig. 3-3 - Fuel Required to Dock and Maximum Coasting Velocity 
vs Initial Range Along -x(l) 
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Fig. 3-4 - Fuel Required to Eliminate Initial Attitude Errors 
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Fig. 3-5 - Fuel Required to Eliminate Initial Attitude Errors 
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Section 4 

CONCLUSIONS AND RECOMMENDATIONS 

A versatile hybrid simulation of orbital maneuvers was developed during 
this study. The program simulates the six-degree-of-freedom motion of a 
chase vehicle as it maneuvers relative to an orbiting target. A cylindrical 
representation of the target is displayed on a CRT display scope so that the 
maneuvers can be visually monitored by a pilot located in the active vehicle. 

Both automatic and manual control can be simulated. A third option includes 
real time switching between automatic and manual control at the pilot's discre- 
tion. System time delays {presently 0.14 second) are included to represent 
delays between commands and thruster firings. Another option permits simu- 
lation of sensor errors so that their effect on accuracy and fuel consumption 
can be assessed. General program application has been maintained where 
possible through modularization and the retention of system parameters as 
program inputs. Computer runs involving pilot interaction must be made in 
real time; however those using only automatic control can be made at accel- 
erated run speeds to achieve a savings in computer time. 

Computer runs using nominal chaser and target data verified the opera- 
tional status of the simulation. Both manually and automatically controlled 
docking runs achieved accuracies within Apollo design limits. Manual attitude 
controllability was improved at the expense of fuel consumption by increasing 
the control torques. Automatic docking generally required much less fuel. A 
two-step docking sequence which yielded low fuel consumption and high accuracy 
was as follows: 

1. Maneuver automatically to a rendezvous point near the docked 
position. 

2. Switch to manual control for final alignment and translation. 


4- 1 


LOCKHEED -HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC D225642-I 


Automatic runs showed the fuel advantage of making maneuvers at low coasting 

« 

velocities. The fuel required to eliminate attitude and attitude rate errors of 
the nominal chase vehicle was also determined. 

Three tasks recommended for continued analysis are described below: 


1. Make docking simulation runs to compare the relative merits 
of purely automatic control, purely manual control, and a com- 
bination of both. Also define desirable switching points between 
automatic and manual control. These runs will require the use 
of hand controllers for manual translation and attitude commands 
to add realism to the manual mode. These controllers need 
supply only on-off capability in the six degrees of freedom. Im- 
provements in the manual mode including minimum impulse bit 
capability and greater control logic sophistication are needed 

to reduce fuel requirements. 

2. Add the capability to simulate remote docking so that free- 
flying unmanned Research and Applications Modules (RAM's) 
can be simulated docking to the space station. This would 
draw heavily from the remote docking work already com- 
pleted during Appendix B of the present contract. 

3. Modify the program to allow simulation of satellite inspection 
and retrieval maneuvers. Include the effects from discrete 
and continuous mass shifts. 
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Appendix A 

HYBRID SIMULATION DIGITAL FLOW CHARTS 
AND SUBROUTINE DESCRIPTIONS 
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Appendix A 


A. 1 GENERAL 

This appendix describes the digital subroutines used in the hybrid sim- 
ulation. Section A. 2 presents a glossary of terms used in this description. 
Variable dimensions are shown in parentheses if they are greater than one. 

The glossary includes numerical values for those variables which are set 
in data statements. 

Sections A. 3 through A.7 describe five of the major digital subroutines 
including subroutine flow charts. In these flow charts I is a free index and 
takes on values of 1, 2 and 3 respectively unless otherwise stated. An 
asterisk (*) preceding a statement identifies it as a hybrid operation. State- 
ments enclosed in quotation marks are messages printed out on the typewriter. 

Simulation timing and interrupts are described in Section A. 8 
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A. 2 GLOSSARY 
Term 

ABRNA, ABRNT 

ADC (16) 

ADSCAL (16) 

ALFAF 

ALT 

ANS 

AST ART 

AX MIS 

BIXYZ(3) 

BRNA, BRNT 

BV(3, 3) 

D(3), DD(3) 

DAC( 16) 

DASCAL(16) 


OF TERMS 


De scription 

jet firing time during automatic portions of the docking 
sequence for attitude control and translation, respectively. 

hybrid variables which are sent from the analog side to 
the digital side. 

scaling factors relating analog and digital values of ADC(l6) 
(see Table A-l). 

scale factor to adjust the display field of view. 

5 

orbital altitude of target (5.0003 x 10 m). 

typed input (yes or no) which determines if the simulation 
run will be automatically stopped at some input time 
(TFLAG). 

number of jet starts during automatic portions of docking 
sequence. 

axial misalignment of the chaser at the end of the docking 
run. 

chaser principal moments of inertia (see Table A-2). 

jet firing time for attitude control and translation, 
respectively. 

transformation matrix from chaser to local vertical 
coordinate s. 

chaser position and velocity relative to the rendezvous 
point expressed in jet axes. 

hybrid variables which are sent from the digital side to 
the analog side. 

scaling factors used in subroutine ALOGIC to relate analog 
and digital values of DAC(l6) (see Table A-l). 


A-2 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC D225642-I 


Term 

DEL(3) 

DISTLV(3) 

DISTB(3) 

DSSCAL(16) 

E(3) 

EDOT(3) 

EPSI(3) 

EYE 

E1DD 

E2DD 

E2S 

F(3) 

FB(3) 

FBRNA, FBRNT 

FCXP, FCYP, 
FSXP, FSYP 

FFA 

FLAG 1 

FSTART 


Description 


maximum change in chaser body rates during one 
digital sample time. 

chaser position relative to the rendezvous point in local 
vertical coordinates. 

chaser position relative to the rendezvous point in body 
coordinates. 

scaling factors used in subroutine HLOOP to relate 
analog and digital values of DAC(16) (see Table A-l). 

chaser attitude errors. 

chaser attitude rate errors. 

Euler angles used to input the initial chaser attitude, 
inclination of target orbit (0 deg). 

attitude error deadband for automatic attitude control. 

used to define automatic attitude control deadzone, 

rate limit during automatic attitude control. 

control forces acting on chase vehicle in local vertical 
coordinates. 

control forces acting on chase vehicle in body coordinates. 

jet firing time for attitude control and translation, respec- 
tively. 

factors used to generate the circle projections in the dock- 
ing target display, 

fraction of TSAMP used for minimum impulse bit jet firing 
a manually controlled flag which terminates the run. 
total number of jet starts (floating point). 
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Term 


Description 


FUEL 

FX, FY, FZ 

FXYZ(3) 

GGCON(3) 

GN(9) 

GSCAL 

HBRNA, HBRNT 

HSTART 

ICHAN 

IDENT 

IDOCK 

IERR 

IFL(6) 

IFL6I 

IGN 

IMAN 


total fuel consumed during docking sequence. 

translation commands (+1, -1, 0) in chaser body co- 
ordinates. These come from the pilot through the 
hand controller during manual control. 

control force per jet along the chaser body axes. 

constants used to calculate gravity gradient torques 
on the chaser. 

current set of gaussian variables from table look-up. 

chaser position scaling factor to permit higher accuracy 
at close range. 

jet firing times during manual portions of the docking 
sequence for attitude control and translation, respectively. 

number of jet starts during manual portions of docking 
sequence. 

channel number (ADC or DAC) which caused flag IERR 
to be set. 

typed input for run identification. 

^+1 successful dock, 
program control flag \ 0 not yet docked. 

\ -1 missed attempt. 

a flag to indicate an overload in the ADC or DAC channels. 

input program control flags (described in Table A-3). 

! 0 = automatic control. 

1 = manual control. 

counter (1-1000) for gaussian variables. 

flag which detects a change in control mode (automatic 
or manual). 
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Term 

INITAL 

IPGMGR, 
JET (I) 

/ JETPR V 

JJT(3) 
MAN CL 
N 

ONE(3) 

OPER 

ORL(3) 

P(3) 

PI 

PMI 

PM 2 

PQRCMD 

QBR(4) 

QIB(4) 


LMSC-HREC D225642-I 


Description 


program control flag (negative the first time through 
subroutine HLOOP). 

ISTEER program control flags. 

state of attitude control jets commanded by automatic 
control laws (+1, 0, -1). 

non-FORTRAN variable used to switch between automatic 
and manual control and controlled by the pilot. 

state of attitude control jets sent to the analog (+1, 0, -1). 

flag which is non-zero during manual attitude maneuvers. 

counter for digital clock. 

upper switching line for manual attitude hold control laws, 
orbital period of target in seconds, 
attitude rate limit under manual control, 
initial chaser body rates (RATEB(3)). 

3.14159. 

range from chaser to target docking collar at which run 
terminates. 

location of rendezvous point along target docking axis 
measured from the target center of mass. 

a flag which is greater than zero whenever attitude changes 
are commanded under manual control. 

quaternions relating chaser attitude to chaser reference 
attitude and used to compute attitude errors during auto- 
matic control. 

quaternions relating inertial coordinates to chaser body 
coordinates. 
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Term 


Description 


QIBD(4) 

desired values of QIB(4) for periods of attitude hold 
under manual control. 

QIV(4) 

quaternions relating inertial coordinates to local vertical 
coordinates. 

QVB(4) 

quaternions relating local vertical coordinates to chaser 
body coordinates. 

QVR(4) 

constant quaternions which describe the target attitude 
relative to inertial or local vertical coordinates. 

RADIUS 

radius of the end circles for the cylindrical representation 
of the docking target. 

RANG, RANGR 

chaser to target range and range rate, respectively. 

RATEB(3) 

angular rates about chaser body axes. 

RE 

£ 

radius of the earth- spherical equivalent (6.36743 x 10 m). 

REF(3) 

Euler angles relating inertial or local vertical coordinates 
to target coordinates. 

RERR, RRERR 

per unit error in measuring chaser to target range and 
range rate, respectively. 

RF(3, 3) 

transformation matrix from inertial or local vertical co- 
ordinates to target coordinates. 

RHOZM, RHOZP 

distances from target center of mass to the rear and front 
faces of the cylindrical representation, respectively. 

RNATT 

number of jets per axis used to produce chaser control 
torques. 

RNTRAN 

number of jets per axis used to produce chaser control 
forces. 

ROLMIS 

chaser roll misalignment at the end of the docking run. 

RPD 

radians per degree (0.017453292). 
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Term 

RSLOPE 

S(3) 

SIGMA(IO) 
SIGR, SIGRR 

SPECI 

STARTS 

TAB(IOOO) 

TAUTO 

TB2T(3, 3) 

TCMAG(3) 

T DELAY 

TFLAG 

TG(3) 

THRUST 

TIME 

TIMEF 

TIMEFI 

TMAN 

TMARG 

TMIN 

TORQN(3) 


Description 

display parameter to size the rear face of the target. 

slope of the manual attitude control law switching curve. 

standard deviations to describe sensor errors. 

standard deviation for range and range rate measure- 
ments, respectively. 

fuel specific impulse (352 sec). 

total number of jet starts (integer). 

input table of gaussian numbers. 

time spent under automatic control. 

transformation matrix from chaser to target coordinates, 
magnitude of chaser control torques. 

system time delay between control commands and firing 
of the jets. 

time at which the simulation run is to be flagged off. 

gravity gradient torques about chaser body axes. 

thrust per control jet (7340N). 

time elapsed since start of docking run. 

desired time to complete docking; must be > TMIN. 

input TIMEF. 

time spent under manual control, 
safety margin on TMIN. 

minimum time required to achieve automatic docking. 

control torque per jet about chaser body axes (see Table 
A-2). 
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Term 

TSAMP 

TWO(3) 

U(3) 

UDOT 
VB(3, 3) 

VMU 

V1M(3), V2M(3) 
V20(3) 

X(3), XDOT(3) 

XB(3), XBD(3) 

XBM(3), XBDM(3) 

XCLZP, YCLZP 

XDOCA, XDDOCA 

XIRL(3) 

XL AT OF 
XLATVE 
XM 


Description 
digital sample time. 

lower switching line for manual attitude hold control 
laws. 

state of the translation control jets (+1, 0, -1). 
target orbital rate. 

transformation matrix from local vertical to chaser 
coordinates. 

gravitational constant times the mass of the earth 
(3.991 x 10 N-m /kg). 

translation dead-zone position and velocity boundary, 
respectively, in each local vertical axis. 

initial chaser velocity (XDOT(3)) in local vertical co- 
ordinates. 

chaser position and velocity, respectively, in local 
vertical coordinates. 

chaser position and velocity with respect to the target 
center of mass in chaser body coordinates. 

measured values of XB(3) and XBD(3) with sensor 
errors included. 

display parameters which locate the front face of the 
docking target. 

chaser position and velocity relative to the target co- 
ordinates. 

inner rate limit for manual attitude control laws, 
chaser lateral offset from the docking axis, 
chaser velocity perpendicular to the docking axis. 

5 

chaser mass (1.26443 x 10 kg). 
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Term 

XQIB 

XSCR(3) 

XSLOPE, YSLOPE 
XSPOT, YSPOT 

Z 


Description 

square root of the sum of the squares of QIB(4); should 
be equal to 1 for orthogonality. 

position of pilot observation point in chaser body coordinate 

display parameters which determine the target orientation. 

display parameters to locate the bright spot on the target 
front face. 

target orbital position at the end of the docking run. 
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A. 3 SUBROUTINE MAIN 

MAIN is primarily a bookkeeping routine. It performs all normal (non- 
error) input and output operations and does much of the initialization required 
before each run. Upon run termination it processes the output data, prints it 
on the in-line printer, and then returns program control to the typewriter. 
Figure A-l shows the flow of subroutine MAIN. Other subroutines which are 
called from MAIN are briefly described below in the order that they are called. 


COEF 

— Computes constant parameters for the automatic 
attitude control logic. 

BEGIN 

— Performs program initialization including the 
computation of initial quaternions and trans- 
formation matrices. 

AUTOIN(D, DD) 

— Initializes automatic translation control logic 
(see Section A.4). 

ALOGIC 

— Sends chaser initial conditions (position, velocity, 
attitude, and attitude rate) to the analog computer. 

HLOOP 

— Performs all digital operations required during 
the simulation, (see Section A. 5) 

MAYMUL 
(O, A, B, C, N) 

— A matrix multiplication routine which computes 
the matrix [C] as follows: 


[c] = [A] [b] 

N is the number of columns in QcJ 


A. 4 SUBROUTINE AUTOIN 

AUTOIN is entered everytime that the simulation is switched into the 
automatic control mode. This routine initializes the automatic translation 
control laws so that docking is achieved at the time originally commanded 
(TIMEFI) if possible. If the time available is too small it attempts to dock 
in minimum time. A flow chart of AUTOIN is shown in Fig. A -2. 
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A. 5 SUBROUTINE HLOOP 

This is the main hybrid subroutine. It completes the initialization of 
each run and then controls all hybrid and digital operations required during 
the simulation. Figure A-3 is a flow chart of subroutine HLOOP. Other 
digital subroutines which are called from HLOOP are briefly described 
below in the order that they are called. 


HLQUAT 

DISPLV 

AUTO 

MANUAL 

TESTSL 
(1, PQRCMD) 

TMFCMD 

AUTOIN(D, DD) 
OUT( 1625) 


computes range, range rate, rendezvous point 
location, and all transformation matrices 
required during the run. 

computes the 11 variables required to generate 
the CRT display. The display equations are 
documented in Ref. 3. This routine also checks 
to see if docking is completed or if a missed 
attempt has occurred. 

computes jet commands under automatic control 
(see Section A. 6). 

computes jet commands under manual control 
(see Section A.7). 

tests for manual attitude rate commands from the 
pilot. 

tests for the control mode (manual or automatic) 
and then tests for manual translation commands 
from the pilot. 

initializes automatic translation control logic 
(see Section A. 4) 

prints out the contents of the flag register. 


The timing of data transfer between the digital and analog computers is con- 
trolled by interrupts in such a way as to produce a system time delay of 0.14 
second. The timing and interrupts are described in Section A.8. 
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A. 6 SUBROUTINE AUTO 

This subroutine computes the commanded control forces and torques 
during periods of automatic control. Figure A-4 is a flow chart of subroutine 
AUTO. The three subroutines called from AUTO are briefly described below. 


ROT LAW (I) 


TJETS(A, B) 


MAYMUL 
(O, A, B, C, N) 


computes the desired state of the torque jets to 
maintain attitude hold. This subroutine was de- 
veloped in a previous study and is fully described 
in Ref. 3. 

computes the desired state of the translation jets 
to achieve docking using the baseline two-burn 
logic. This subroutine was developed in a pre- 
vious study and is fully documented in Ref. 2. 

matrix multiplication routine described in Section 
A. 3. 


The torque commands (JJT(I)) and force commands (F(I)) are delayed by 0.14 
second before being sent to the analog computer. 


A. 7 SUBROUTINE MANUAL 


MANUAL computes the torque commands (JJT(I)) during periods of man 
ual control according to the switching logic described in Section 2.5. A flow 
chart of this subroutine is presented in Fig. A- 5. 

A. 8 SIMULATION TIMING AND INTERRUPTS 


During the simulation program variables are transferred between the 
digital and analog computers in the following ways. 

1. Sixteen digital-to-analog (DAC) conversion channels. 

2. Sixteen analog-to -digital (ADC) conversion channels. 

3. Sixteen logical data bits (+ or 0) from digital to analog through the 
data word. 

4. Sixteen logical data bits (+ or 0) from analog to digital through the 
data word. 

5. Eight hard-wired sense lines. 
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This data transfer is pictured in Fig. A-6. Three digital subroutines involved 
in the transfer of data (TMFCMD, JETFIR, and TESTSL) are also pictured 
and described below. 

TMFCMD interrogates the data word from the analog for control 

mode (manual or automatic) and manual translation 
commands. 

JETFIR sends the attitude jet firing commands to the analog via 

the data word; also processes all jet firings to calculate 
the number of jet starts, translation burn periods, and 
rotational burn periods. 

TESTSL interrogates the sense lines for manual attitude rate 

commands. 

Simulation timing is controlled from the digital by the two interrupts 
shown in Fig. A-7. The timing sequence for the first 0.54 second of the sim- 
ulation is shown in Fig. A-8. Clock interrupts F. and T. represent interrupts 

"FROM" and "TO," respectively. 

During F^ (at time = 0 ) the control status (manual or automatic) and 

the manual force commands (MFCMD) are sampled in subroutine TMFCMD. 
The initial values from the analog are stored in the ADCs at this time also. 
The ADCs are inhibited from storing new analog values until the ADC-EBL 
line (on the analog) goes high. This occurs every 100 milliseconds (after 
t=0). The ADC-EBL pulses are obtained by counting 100 pulses from the 
analog 1 kc clock source. The width of the ADC-EBL pulses is determined 
by the time required to store the analog values in parallel and sequentially 
convert to digital information in the main stream of subroutine HLOOP. The 
clock interrupts initiate sampling of the analog information (via the ADCs and 
data word) and sending information (via DACs and data word) to the analog. 

In order to sample MFCMD at the same time the analog values are 
stored in the ADCs (interface), an initial time of 100 milliseconds is counted 
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first. When interrupt occurs the control status (CS) and MFCMD are 

sampled in subroutine TMFCMD. Then 40 milliseconds is loaded into the 
digital timer. After the 40 milliseconds interval, interrupt T ^ occurs and 

the computed information from F^ is sent to the analog in subroutine JETFIR. 

The effect of this timing is to produce a delay of 140 milliseconds between the 
commands and their effect on the vehicle dynamics. The timer is then loaded 
with 60 milliseconds so that the next interrupt (F^) will occur at the same time 
that the ADC-EBL line goes high to sample new analog information. The inter- 
rupts occur every 40, then 60 milliseconds until the simulation is completed or 
aborted. This coincides with the 100 millisecond ADC-EBL pulses. 

There is no conflict in timing if the digital computation can be completed 
in less than 100 milliseconds because the digital cannot execute the instruction 
to store the analog values until the ADC-EBL line goes high. Therefore, the 
digital will wait until it can execute the store instruction and then proceed with 
the computation. The time to complete the computation is typically less than 
40 milliseconds. A third interrupt which operates in HLOOP is shown below. 
This interrupt is executed anytime a timing error occurs. 


TIMING ERROR 
INTERRUPT 
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Table A-l 


NOMINAL VALUES FOR ANALOG-TO-DIGIT AL 
AND DIGIT AL-TO -ANALOG SCALING 


ADSCAL(I) 


DASCAL(I) 


DSSCAL(I) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 


0.1 

0.1 

0.1 

50.0 

50.0 

50.0 

-5.0 

-5.0 

-5.0 

2.0 

2.0 

2.0 

2.0 

0.1 

0.1 

0.1 


0.00001 
0.000005 
0.000005 
0.025 
0.025 
0.025 
0.00005 
0.00005 
0.00005 
0.0 
0.0 
0.0 
0.0 
- 0.0 
0.0 
0.0 


1.0 

1.0 

1.0 

- 1.0 

0.025 

0.025 

0.00005 

0.00005 

0.00005 

1.0 

1.0 

- 1.0 

1.0 

1.0 

- 1.0 

- 1.0 


Table A-2 


NOMINAL VALUES FOR CHASER PRINCIPAL MOMENTS OF INERTIA 
AND CONTROL TORQUES PER THRUSTER 



BIXYZ(I) 

TORQN(I) 

I 

(kg-m 2 ) 

(N-m) 

1 

4.203 x 10 6 

1.3558 x 10 4 

2 

2.052 x 10 7 

2.8811 x 10 4 

3 

2.319 x 10 7 

2.5421 x 10 4 
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Table A- 3 

PROGRAM OPTIONS BY INPUT FLAGS - IFL(I) 


FLAG 

Option Description 

Input 

Value 

IFL(l) 

Target in local vertical hold 

1 


Target in inertial hold 

2 

IFL (2) 

Perfect sensors 

0 


Effect of sensor inaccuracies included 

1 

IFL (3) 

Control torques = 0 

0 


Baseline automatic attitude control 

1 

IFL (4) 

Control forces = 0 

1 


Translation jets aligned with local vertical axes 

2 


Translation jets aligned with body axes 

3 


Same as above and including sensor errors 

4 

IFL (5) 

Not currently used 


IFL (6) 

Automatic control 

0 


Manual control 

1 
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Fig. A-l - Flow Chart of Subroutine MAIN 
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* SELECT INITIAL CON DITI0N 

* CLEAR ALL DAC REGISTERS 

* SET C0NTR0L LINE 0 

* RESET C0NTR0L LINE 1 

* SET C0NTR0L LINE 6 

* STEERING AND CONSOLE SELECTION 

* ZER0 DATA W0RD T0 ANAL0G 



} 



IFL(6) = IFL6I 



i 



Print: IDENT, IFL, X, V20, EPSI, REF, P, TIMEF 






CALL BEGIN 

No 


Yes 

*' 



17 CALLAUTC 

<IN(DISTLV, XD0T) 

I ;^i 

18 CALL AUT0IN(DISTB, XBD) 


1 


1 



* v ) 



XQIB = SQRT(QIB(1)**2+QIB(2)**2+QIB(3)**2+QIB(4)**2) 






Print: QIB 

, XQIB, QIV, QVB, QBR, QVR 



CALL AL0GIC 
CALL HL00P 



Z = 2.*HU 

XQIB = SQRT(QIB(1)**2+QIB(2)**2+QIB(3)**2+QIB(4)**2) 
FBRNA = BRNA*TSAMP 
FBRNT = BRNT* TSAMP 
FSTART = STARTS 

FUEL = (THR UST/ SPECI)*(RNTRAN*FBRN T+RNATT* FBRNA) 


Yes 

X 

LAN 

\ No 


^ ^ 

TMAN = TIME-TAUT0 


TAUT0 = TIME - TMAN 

ABRNA = ABRNA * TSAMP 


HBRNA = HBRNA * TSAMP 

ABRNT = ABRNT * TSAMP 


HBRNT = HBRNT* TSAMP 

HBRNA = FBRNA -ABRNA 


ABRNA = FBRNA - HBRNA 

HBRNT = FBRNT -ABRNT 


ABRNT = FBRNT -HBRNT 

HSTART = FSTART - ASTART 


ASTART - FSTART - HSTART 

♦ r (°lt * 


Fig. A 

'I > 

(Continued) 
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Fig. A- 1 - (Concluded) 
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Fig. A-2 - Flow Chart of Subroutine AUTOIN (D, DD) 
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©*— C 



M2) 


0 ^ '» (k>oo) 


* CLEAR ALL INTERRUPTS 

* STORE ADDRESS OF "INTRTP2 + 1" IN '62 

* SET ENABLE FLAG 

* ENABLE EXTERNAL INTERRUPTS 

* STOP INTERROGATING INTERRUPTS 
CALL TESTSL (I, PQRCMD) 

INITAL = 0 

* BL0CK TRANSFER DAC'S 0 — 1 5 

* RESET FLAG 1 

* HALT PR0CEDE ON "EXECUTE" 

* START INTERROGATING INTERRUPTS 

* ST0RE ADDRESS OF "FR0M+1" IN '46 

* ST0P DIGITAL TIMER 

* L0AD TIMER WITH 100 MILLISEC 

* START DIGITAL TIMER 
CALL TMFCMD 

* PERMIT EXTERNAL INTERRUPTS 




* SELECT 0PERATE 

* RESET ENABLE FLAG 

* LOAD ADCU6) INTO INTERFACE 

* WAIT FOR TIMING INTERRUPT 
(all normal interrupts occur here) 


3 

IFL(6) = "/JETPRV" 1 


TAUT0 

IMAN 

FSTART 

FBRNA 

FBRNT 

ABRNA 

ABRNT 

ASTART 


TIME - TMAN 
1 

STARTS 

BRNA 

BRNT 

FBRNA - HBRNA 
FBRNT - HBRNT 
FSTART - HSTART 



TMAN = TIME - TA UT0 
IMAN = 0 

FSTART = STARTS 
FBRNA = BRNA 
FBRNT = BRNT 
HBRNA = FBRNA -ABRNA 
HBRNT = FBRNT -ABRNT 
HSTART = FSTART -ASTART 



Yes 


IFL(4)> 2 


» ^~CALL AUT0IN(DISTB, XBDI 


| CALL AUT0IN (DISTLV, XD0T) 


I FX = FY = FZ = 0 I 



1 12lW— 


Fig.A-3 - Flow Chart of Subroutine HLOOP 
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Fig.A-3 - (Continued) 
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CALL HLQUAT 

TG(1 ) = GGC0N(l) * VB(2, 1) * VB(3, 1) 
TG(2) = GGC0N(2) * VB(1. 1)* VB(3, 1) 
TG(3) = GGC0Ni3) * VB(1, 1) * VB(2, 1) 
CALL DISPLV 


RESET ENABLE FLAG 
FB(1 ) = FX*FXYZ(1) 
FB(21 = FY * FXY Z(2) 
FB(3) = FZ* FXYZ(3) 
CALL MANUAL 


~C 


Yes 



CALL AUT0 

FX = U(I)* (ABS(U{1) + .2) ) 
FY = U(2) * (ABS(U(2) + .2) ) 
FZ = U< 3) * ( ABS( U ( 3) + .2) ) 


DAC(l) 
DAC (2) 
DA C ( 3 ) 
DAC(4) 
DA C ( 5 ) 
DAC (6) 
DAC (7) 
DA 08) 
DAC(9) 
DAC(iO) 
DAC(ll) 
DA C ( 1 2 ) 
DA C ( 1 3 ) 
DA C ( 1 4 ) 
DA C ( 1 5 ) 
DA C ( 1 6 ) 


DSSCAL(l) * 
DSSCA L(2) * 
DSSCA L( 3) * 
DSSCA L (4) # 
DSSCA L( 5) * 
DSSCA L (6) * 
DSSCA L(7) * 
DSSCA L(8) * 
DSSCA L(9) * 
DSSCA L( 10) 
DSSCA Lll) 
DSSCA L( 12) 
DSSCA L(1 3) 
DSSCA L(1 4) 
DSSCA L(1 5) 
DSSCA L( 16) 


XSL0PE 

YSL0PE 

XSP0T 

YSP0T 

TG(2) 

TG( 3) 

F(1 ) 

F(2) 

FI 3) 

* RSL0PE 

* FSXP 

* FCXP 

* FSYP 

* FCYP 

* XCLZP 
* YCLZP 



* L0AD DAC(16) INT0 INTERFACE 

* RESET ENABLE FLAG 



* SELECT H0LD 

* DISABLE EXTERNAL INTERRUPTS 
ICHAN = J+3 

Type: ICHAN, TIME. (DAC(.T), .1=5, 9) 
ISTEER = 4 




Accept: IPC'.MOR 



* SELECT 0PERATE 

* ENABLE EXTERNAL INTERRUPTS 


ISTEER = 3 
* SELECT H0LD 
CALL 0UTI162S) 



I * SELECT H0LD 
f RETURN J 


Fig. A-3 - (Cone luck'd) 
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SUBROUTINE AUTO 


"TFL(2)> (T 


SIGR = SIGMA (5) * RANG 


SIGR < SIGMA (4? 


SIGR = SIGMA (4) 


SIGRR = SIGMA(7) * RANGR 


SIGRR < SIGMA ((Tf 


SIGRR = SIGMA (6) 


r D0 77 1 = 1,9 

I IGN = IGN + 1 


S IGN>J1009 


! I GN(I) = TAB(IGN) 

77 C0NTINUE 


IFL<3) > 0 


ED0T(1 ) = RATEB(1 ) + SIGMA(8) * GN(6) 

ED0T (2 ) = RATEB(2) + SIGMA{9)*GN(7) 

ED0T(3) = RATEB{3) + SIGMA (10) * GN<8) - UD0T 
-|- D0 70 1 = 1,3 

I E(I) = -2. *QBR(I+1) -SIGMA(I)*GN(I) + TDELAY*ED0T(I) 
, CALL R0TLAW(I) 

I JJT(I) = JET(I) 

70 C0NTINUE 




CALL TJETS(DISTLV, XD0T) 
F(I) = U(I) * FXYZ(I) 

CALL MAYMULfO, VB, F, FB, I) 


RERR = SIGR * GN(4) CALL TJETS(DISTB, XBD) 

RRERR = SIGRR *GN(5) FB{I) = U(I)*FXYZ(I) 

XBM(I) = DISTB(I) + RERR * DISTB{I)/RANG CALL MAYMULfO, B V, FB, F, 1 ) 
XBDM(I) = XBD(I) + RRERR*XBD(I)/RANGR 
CALL TJETS(XBM, XBDM) 

FB(I) = U(I) * FXYZ(I) 

CALL MAYMULfO, B V, FB, F, I) 
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Fig. A-5 - Flow Chart of Subroutine Manual 
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RESET ENABLE FLAG 

STOP DIGITAL TIMER 

LOAD TIMER WITH N 

START DIGITAL TIMER 

STORE FLAG REGISTER IN 3003S 

STORE ACCUMULATOR IN 3050S 

STORE ADDRESS OF "TO + 1" IN '46 INTERRUPT 

"FROM" 

CALL TMFCMD 
N = 6000 (60 millisec) 

RESTORE FLAG REGISTER FROM 3003S 
RESTORE ACCUMULATOR FROM 3050S 
RESET INTERRUPT 
ENABLE INTERRUPTS 
RETURN TO "MAIN STREAM" 



RESET ENABLE FLAG 

STOP DIGITAL TIMER 

LOAD TIMER WITH N 

START DIGITAL TIMER 

STORE FLAG REGISTER IN 302 IS 

STORE ACCUMULATOR IN 3050S INTERRUPT 

"TO" 

CALL JETFIR 

BLOCK TRANSFER DAC'S 0— *-15 
N = 4000 (40 millisec) 

STORE ADDRESS OF "FROM+1" IN '46 
RESTORE ACCUMULATOR FROM 3050S 
RESTORE FLAG REGISTER FROM 302 IS 
RESET INTERRUPT 
ENABLE INTERRUPTS 
RETURN TO "MAIN STREAM" 


Fig. A-7 - Normal Timing Interrupts 
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l| II 

ADC-EBL Pulses 

n fi fi 

II 

n n 

H 

0 0 

U 

__JL_ 

_JL 

F FT 

1 2 1 

Clock Interrupts 

F 3 T 2 

1 1 

F 4 T 3 

F 5 T 4 

| 1 

F 6 

1 

T 5 

1 

| 1 ^ 

0 .1 .14 

1 I 

.2 .24 

.3 .34 

.4 .44 

.5 

.54 


Time (sec) 


Fig. A-8 - Timing Sequence for the First 0.54 Seconds 
of Hybrid Simulation 
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Appendix B 

ANALOG EQUATIONS 
AND WIRING DIAGRAMS 
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Appendix B 

B.l ANALOG EQUATIONS 

The chaser equations of motion (translation, rotation, and kinematic) 
are shown below as they are solved on the analog computer. In these equa- 
tions a dot over a variable indicates its time derivative (i.e., x = position, 
x = velocity, x = acceleration). A quantity in brackets, [] , is a scaled 
analog term where 1.0 equals 100 volts. 


Translation 


[±m 

[0.5 

= +200 00 

x (2) 

+ 6000 0) 2 

x (1) 

50 

1000 

[ x (2) 

= -20oo 

i(l) 

, 40,000 


[ 0.5 

5 

+ M 

20,000 


40,000 

M 


20,000 


x (3) 
0.5 


-2000O) 2 


x (3) 

. 40,000 

F(3) 

1000 

+ M 

20,000 


where, 

x (i) = chaser position in local vertical coordinates (m) 

M = chaser mass (kg) 

0) = target orbital rate (rad/ sec) 

F (i) = control forces expressed in local vertical coordinates (N) 
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Rotation 


Automatic Mode 



Manual or Switching Mode 



[lOE(l) 

- 

10P 


1 

!lp + 

j 


lOPcj 

♦1 

LP-! 


iopc] 

[lOE(2) 

= 

ioq] 

- ' 

|lq+] 

! 

1 

riooc] 

+ 1 

} L Q-j 


^lOQcj 

[lOE(3) 

= 

lH 

1 _ 1 

LR +j 


[lO Rcl 

+ ' 

| 

1 TO _j 
|L R | 

. 

|^10RcJ 
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where 

P.Q.R 

PC, QC, RC 

^xx’ ^yy’ ^zz 
TC(i) 

TG(i) 

T , T , T = 

x y z 

LTC1+, LTC2+, 
LTC3 + 

LP+, LQ+, LR+ = 
E(i) 


angular rates about chaser body axes (rad/sec) 

(same as RATEB(i) defined in Appendix A) 

angular rate commands from pilot 

chaser principal moments of inertia (kg - m ) 

control torques about chaser body axes (N - m) 
gravity gradient torques about chaser body axes (N - m) 
absolute value of TC(i) 

logical variables defined in Table B-l 
logical variables defined in Table B-2 

chaser attitude rate errors (rad/sec). 


Illegal combinations of logical variables are prevented from occurring. 
Illegal sets are those which require opposing jets to fire and result in no n6t 
thrust (i.e., LTC1 + = LTC1 - = 1, LP + = LP - = 1, etc). 


Kinematics 


[lO QIB(l) | = - 1 5 P QIB(2)J - |5Q QIB(3) j - j^5 R QIB(4) 
10 QIB(2) J = + |5P QIB(l)j - ^5Q QIB(4) j + |5R QIB(3) 


10 QIB(3) I = + 


5 P QIB(4) 


1 


5 Q QIB(l) 


10 QIB(4) 


5 P QIB(3) + 5Q QIB(2) 


QIB(i) 


2 V QIB(1) 2 + QIB(2) 2 + QIB(3) 2 + QIB(4) 


1 

|5R QIB(2) J 
+ [5R QIB(l) J 

[QXB(i)] i=1 , 2 ,3. 


where, 


QIB(i) = quaternions relating inertial coordinates to chaser body 
coordinates 

QIB(i) = orthogonalized QIB(4) 
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Table B- 1 

LOGICAL VARIABLES FROM DATA WORD WHICH 
REPRESENT ATTITUDE THRUSTER FIRINGS 


LTC1+ = 1,0 (logical) for positive attitude thrust (P channel) 

LTC 1 - = 1,0 


negative 


(P channel) 

LTC2+ = 1,0 


positive 


(Q channel) 

LTC2 - = 1,0 


negative 


(Q channel) 

LTC3+ = 1,0 , 


positive 


(R channel) 

LTC 3 - = 1,0 (logi 

cal) for negative attitude thrust (R channel) 


Table B-2 

LOGICAL VARIABLES WHICH REPRESENT RATE 
COMMANDS FROM THE PILOT 


LP+ = 1,0 (logical) 
LP - = 1,0 
LQ+ = 1,0 
LQ - = 1,0 
LR+ = 1,0 , 

LR - =1,0 (logical) 


for positive rate command 


negative 


positive 


negative 


, positive 



for negative rate command 


(P channel) 
(P channel) 
(Q channel) 
(Q channel) 
(R channel) 
(R channel) 
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B.2 ANALOG WIRING DIAGRAMS 

Symbols used in the analog wiring diagrams are defined in Fig. B-l. 
Diagrams of the translation, rotation, and kinematic equations are shown in 
Figs. B-2, B-3, and B-4, respectively. 

In Fig. B-2 the analog switches SW001, SW011, SW021, and SW031 are 
switched left (L), facing the panel, to decouple the translation equations and 
they are switched right (R) to add cross coupling. Control line 1 (CL1) is used 
to scale up the analog position values for improved accuracy when all three 
axes are less than 45 meters. 

The manual command logic is shown in Fig. B-5; display logic is shown 
in Fig. B-6, and timing logic is shown in Fig. B-7. 

B.3 AUTOMATIC ANALOG SET-UP 

A hybrid program is used during set-up of the simulation. This pro- 
gram calculates the analog potentiometer values and automatically sets the 
potentiometers. It also provides the operator with a check list of required 
setup procedures. A listing of this program is presented in Fig. B-8. 
Standard symbols are used where possible. Other symbols not previously 
defined are: 

W1 = CO 

W2 = CO 2 

J1 = label of a fixed value 

„ C002 = address of potentiometer No. 002. 
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Summer 





-O- 

Potentiometer 



* 



(1*B + 10* A)dt 


D = - (1* A+ 10* B) 


D = ^* (1* A+10*B) 


D = -A if T = 1 (logical) 
D = -B if T = 0 (logical) 



Digital Control Line 6 


Fig. B-l - Symbols Used in Wiring Diagrams 
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Digital-to -Analog Converter Output 



Analog -to -Digital Converter Input 


1 



Normally Open Contact 

Normally Closed Contact 
Logical Relay Control 




D = - 1 .0 * B if L = 0 (logical) 

D = -1.0* B - 10* A if L = 1 (logical) 



D = A • B 
E = A • B 


Fig. B-l - (Continued) 
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Logic Trunk 


L = 1 if A > 0 volts 

L = 0 if A < 0 volts 


Comparator 


D = 1 if (A + B) > 0 volts 
D = 0 if (A + B) < 0 volts 




Resolver 


Fig.B-1 - (Concluded) 

j 
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©- 
T4 b\ 
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Fig. B-2 - Analog Wiring Diagram for Chaser Translation 








Fig, B-3 - Analog Wiring Diagram for Chaser Rotation 
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Fig. B-3 - Concluded 
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Fig. B-4 - Analog Wiring Diagram for Chaser Attitude 
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Fig. B-7 -Analog Wiring Diagram of Timing Logic 
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I JOB H-10 HEESCHEN N8-5C99-1003 2/1/72 

5 H 0 1 , , I U 1 * 

JCE.XMODE? aC0NTINUE-0N-ERR0R OPTIONa 
1 , CONSOLE ) 

JPS.MODEl 

aSPACE SHUTTLE DOCKING STUDYaJ! aMKSaI 
aSET DIGITAL T! HER J REDsREDa* 
aSET ANALOG TIMER! A*22, B*22» C*0S, D* 08 a J 
aSET DFG’S TO INV.aJ 

I NPUT > a ARE DAC»S OUT OF LOOP TEST?a*aTYPE GOaI 
1,0100) AlNPUT DATA a 

1.0110) BC=. 0090, B0=, 0000. 

1,0200) Wl=l, 1067736-3, W2*l , 2249465E-6 , 

1,0300) M = 1 , 264426E+5 , 

1,0310) TX » 2 , 71166+4 , 

1.0320) TY * 5.76226+4, 

1,0330) TZ * 5,08426+4. 

1.0400) IXX * 4 , 2030E+6 , 

1,0500) IYY « 2 , 0520E+7 , 

1.0600) IZZ * 2,31906+7, 

1,0610) PC=. 3500, QC=, 3500, RC*. 3500, 

1.0700) Jls.0100, J2=, 1000, J3*. 2500, J4s, 4000, J5e, 0500 jj6«, 8300, 

1,0710) J7 = ,4400, J8=. 2380, J9=. 4570, J10s, 2000, 

1.0720) aCOEFFICIENT CALCULATJONa 

1 . 0800 ) C 000 = J5,C001a6.E+4#W2, 0003=40000, /M,C010aB0,C011=BC# 

1,0810) C0Pl2a(100.*TY)/lYY,C030aj7,C03l*<10.*TX)/IXX,C032a<10.*TX)/IXX, 

1,0900) C020sl000.*W2,C02l a J10;C022aB0,C023.BC;C101aJ2; 

1,0910) Cll0a(l00.*TY)/IYY,Cll?a<l00,#TZ)/IZZ,Cl22*(l00,* T z>/izz, 

1,1000) 0200=40000, /M, C201aj5,C202aB0,C203a200,*W*,C2l0aj5,C211ag2, 

1,1100) 0212=4, E+5/IXX, 0222*80, C 25 3»J3, 03 00 =20, *Wl,C302a2,E+6/IYY, 

1,1110) C230ajl,C240sjl,C250aJl, 0012=40000 , /M , 

1,1120) C050»PC,C220*RC, 0231*00, 

1,1200) C303sJl,C3l0*(IZZ + 'IYY)/IXX,C3ll*Jl,C3l2*l,E + 6/IXX,C322*BC, 

1,1300) 0351 = J 4, 0400 = J2, 0401 8 B0,C402b 4,E + 5/ IYY, 0410 = J2iC4^.i»B0^ 

1,1400) C412=4.E+5/IZZ,C413aBC,C500aBC,C501=( IZZ-lXX)/IYY,C5l0=Ji; 

1,1500) C511 = 2.E + 6/IZZ,C5l2a< I YY- 1 XX ) / 1 ZZ , 0600* J2 , 0601* J2 , 0610* J2 , 

1,1600) C611«J2,C700aJ8,C040Bj9,C041*J9, 

JDWL, OUTPUT J 
1? AL I ST PROGRAMa 

lj aEXECUTE PROGRAMa 

tCOEF : aLIST COEFFICIENTSa 

JS,/T, OUTPUT) 

JPS,MODE?$COEF,SET| a ALL POTS SETa* 

JS,JT, INPUT? aPAUSE TO CORRECT POTSa*aTYPE GOa? 

,0004, VERIFY JICOEF.READI aALL POTS CHECKEDaI 
NORMAL ? 

JCL00 , 0 , SET ? 

JCL01,0,SETJ 
701,06,1, SET? 

JS, JT, INPUT?aTYPE GOa* 

3, CONSOLE? aEX I T PROGRAMa 

$ JOB 

SRW , ,SU1 

Fig.B-8 - Listing of Hybrid Set-up Routine 
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